Effect of exchange interaction on superparamagnetic relaxation 
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We use Langer's approach to calculate the reaction rate of a system of two (classical) spins 
interacting via the exchange coupling J in a magnetic field H, with uniaxial anisotropy of constant 
K. We find a particular value of the exchange coupling, that is j = J/K = j c = 1 — ft 2 , where 
h = H/2K, which separates two regimes corresponding to a two-stage and one-stage switching. For 
j jc the Neel-Brown result for the one-spin problem is recovered. 



PACS numbers: 05.40.-a, 75.10.Tt 



I. INTRODUCTION 



Due to their high coercivity, single-domain nanoparticles exhibit long-range stability of the magnetization and 
thereby that of the information stored in recording media. The storage density may be increased by using very small 
particles, but then surface effects become dominant and affect the magnetization relaxation. Therefore, one of our 
aims is to include surface effects in the calculation of the relaxation time of the particle magnetization. However, this 
requires a microscopic approach to account for the local environment inside the particle, and thus include microscopic 
1 j interactions such as spin-spin exchange, in addition to anisotropy and applied field. Unfortunately, this leads to 
' O ■ a rather difficult task owing to the large number of degrees of freedom which hinders any attempt to analyze the 
energyscape. For this reason, inter alia, calculations of the reversal time of the magnetization of fine single-domain 
ferromagnetic particles, initiated by Neel 0, and set firmly in the context of the theory of stochastic processes by 
Brown |2j, have invariably ignored all kind of interactions and included only the internal anisotropy of the particle, 
the random field due to thermal fluctuations, and the Zeeman term. 

In 1968 Langer jj] came up with a general approach to the calculation of the reaction rate of a system with N 
degrees of freedom, which is a powerful method of attack on the problem of interaction. Within this approach the 
problem of calculating the relaxation time for a multi-dimensional process is reduced to solving a steady-state Fokker- 
Planck equation (SSFPE) in the immediate neighborhood of the saddle point that the system crosses as it goes from 
a metastable state to another state of greater stability. The idea is that a steady-state situation can be set up by 
7-H ' continuously replenishing the metastable state at a rate equal to the rate at which it is leaking across the activation 
energy barrier. A brief review of the main steps and application to the problem at hand are given after Eq. JSJ), see 
also 4] for uniform and non-uniform magnetization, and 5] for comparison with Kramers' theory. However, Langer's 
approach is only valid in the limit of intermediate-to-high damping (IHD) because of the inherent assumption that 
the potential energy in the vicinity of the saddle point may be approximated by its second-order Taylor expansion. 
The result for small damping fails because the region of deviation from the Maxwell-Boltzmann distribution set up 
. in the well extends far beyond the narrow region at the top of the barrier. 

The problem associated with the generalization of Brown's theory to include interactions is rather difficult and 
can in general be solved only numerically. But before attacking this problem, one needs to understand the effect 
O ■ of exchange interaction on the relaxation time of the minimal system, i.e., for a pair of atomic spins coupled via 
exchange interaction, including of course the usual magneto-crystalline anisotropy and Zeeman terms. Besides, in the 
case of exchange interaction, this is the unique non-trivial step towards the above-mentioned generalization, where 
analytical expressions can be obtained for the relaxation time. It is the purpose of this work to solve this problem 



within Langer's approach and to compare with the Neel-Brown result for the one-spin problem. As a byproduct, 
we will show that Langer's quadratic approximation at the saddle point fails when the exchange coupling assumes a 
"critical value" even in the IHD limit. In this case, finiteness of the relaxation rate requires a higher-order expansion 
of the energy near the saddle point. The generalization to a multi-spin particle will be briefly discussed at the end. 
We finally mention that the effect of dipolar interactions on the relaxation time has been considered by a few 
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FIG. 1: Energyscape from Eq. iJTJ (tp = 0) for j = 0.25, 0.8, 1.5 and h = 0.1. The arrows indicate the switching paths. 

authors. In Refs. [|| approximate analytical expressions were obtained for the relaxation time in various situations of 
assemblies of magnetic moments. In Ref. @ a pair of coupled dipoles was dealt with using the (numerical) Langevin 
approach. 

II. STATEMENT OF THE PROBLEM AND ANALYSIS OF THE ENERGYSCAPE 



We consider a system of two exchange-coupled spins with the Hamiltonian 

(Si • ei) 2 + (S 2 • e 2 ) 2 l - h ■ (Si + S 2 ), 



n = - ] -s 1 -s 2 - l - [ 



j = J/K, h = H/2K 



(1) 



where J > is the exchange coupling, K > the anisotropy constant, H the applied magnetic field, and h the reduced 
field, i.e., < h < 1. ffj are uniaxial anisotropy unit vectors. Here we restrict ourselves to the case h \\ ei,i = 1,2. 
Owing to the symmetry of this system with respect to rotations around the easy axis, the number of variables reduces 
to three, 9\ , 9 2 , <p = (fi — <^2 ■ 

Now we apply Langer's approach to the energy Q and study the relaxation rate as a function of the exchange 
coupling j. We first analyze the energyscape in Fig. ^ The absolute minimum of the energy JQ) corresponds to the 
ferromagnetic order along the easy axis, 



(0,0, p); e 



(0) 



J 



-2/i-l, 



(2) 



where henceforth denotes the energy of the state. One metastable state corresponds to a ferromagnetic order 
opposite to the field 



There is also the metastable state of anti-ferromagnetic order, 



J 



2/i-l. 



(0,7r ;¥ p)or(7r,0; V ); e<°> = | - 1. 



(3) 



(4) 



As to saddle points, we find that their number and loci crucially depend on the exchange coupling constant j. More 
precisely, for j > j c = 1 — h?, there is a single saddle point given by 



(cos6>i = cos6> 2 = -h; ip); e^ 0) =-| + /i 2 , 
whereas for j < j c there are two saddle points given by 

cos ^,2 = \ (~h-a £ ± V&) =Xg; if, 



(5) 



(6) 
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where e = ± and a £ = ey/1 - j, b 6 = e-y/l + (j/2) 2 /(l - j), A e = (h + a £ ) 2 + 2j - Wh, with energy, 



(°) = -L / (1 -4- fce/>. - ^2 _ C„e r M2 _ jfi 4- fc £ /). _ + If 



2 y (1 + b^h ±y - (a* + hf - ±(l + Vh -±) + -(h 2 - (a^f + 2b*h). (7) 

At j — j c the saddle point 10 with e = + merges with the saddle point J5J , while that with e = — merges with the 
metastable state (@J , see Fig. ^ central panel. 

Starting with both spins aligned in the metastable state 10 with #1.2 = 7T, if j < j c one of the two spins crosses the 
saddle point (©(£ = +) into the state by reversing its direction. Then the second spin follows through the second 
saddle point (JSJ(er = — ), of lower energy due to the exchange coupling (see Fig. (2J (left)), ending up in the stable state 
@. In Fig. Q] (left) the path is indicated by a pair of curved arrows. There are actually two such paths corresponding 
to the two-fold symmetry of the problem owing to the full identity of the two spins. Note that when the first spin 
starts to switch and arrives at 0\ ~ 7r/2, the second spin has #2 ;$ tt (hence the curved arrows in Fig. [IJ, which 
suggests that in the switching process of the first spin, the position of the second spin undergoes some fluctuations 
creating a small transverse field, and when B\ = the second spin goes back to the position 62 = tt before it proceeds 
to switch in turn. The successive switching of the two spins through the corresponding saddle points is a sequential 
two-step process [T^ |. so the relaxation rates for j < j c add up inverse- wise. In the case j > j c the two spins cross the 
unique saddle point (JSJ to go from the metastable state © into the stable one J2J in a single step, see Fig. 2] (right) 
where the path is indicated by a single straight arrow. There is the symmetry ±#( s ), which leads to a factor of 2 in 
the relaxation rate. Therefore, if we denote by F^^. , Tj K j , and Tj>j c the respective relaxation rates, the relaxation 
rate of the two-spin system is given by 

r = J 2r /<ieW( r ;<; c + r 7<;J> *i<io (8) 

l2Ti>i.. if j>j c . 

In Langer's approach the Ts are obtained from the SSFPE linearized around each saddle point, using the fluctuating 
variables r\i = (U,Pi) where tj = 0, — 9f, pi = (fi — ip*, i = 1,2, or more adequately the "canonical" variables 
4>n = (C±jC±) w ith C± = (ii±*2)/v2j C± — (Pi =tP2)/\/2, in which the energy Hessian is diagonal. The deterministic 
dynamics of the system is governed by the Landau-Lifshitz equations, which upon linearization near the saddle point, 
read 

d t U = -d Pz TL 2 - adtiV.2 ^ 
d tPi = -ad Pi H2 + d ti H 2 , i = l,2 [y> 

where TI2 is the quadratic approximation of the energy Q near the saddle point and a is the damping parameter. 
In matrix form Eqs. © become dtrji — Y) j Mij d rij Ti.2, where M is the dynamic matrix containing the precessional 
and dissipative parts. Rather than investigating the stochastic trajectories ?7i(t) that arise by adding a noise term 
in Eqs. ©, Langer concentrates on the distribution function p({r]},t) as the probability that the system is found 
in the configuration {77} at time t. The time evolution of p is governed by the Fokker-Planck equation (FPE) 
dtp + ^2id ni Ji = 0, which is a continuity equation with the probability current Ji given, near the saddle point, by 
Ji = — Mij (d rjj Ti.2 + (l/fceT)^.) p. In order to calculate the nucleation rate, one must solve the FPE. A particular 
solution is obtained at equilibrium and is given by the Maxwell-Boltzmann distribution p eq — exp(— /37i {rj})/Zo, which 
corresponds to zero current, Ji = 0. However, what is really needed is a finite probability current flowing across the 
saddle point. In fact, instead of solving the time-dependent FPE, Langer solves the SSFPE dtp — near the saddle 
point. The steady-state situation may be realized by imposing the boundary conditions: p ~ p eq near the metastable 
state and p ~ beyond the saddle point. Then, the problem of calculating the escape rate reduces to the calculation 
of the total current by integrating the probability current Ji over a surface through the saddle point. This led Langer 
to his famous expression for the relaxation rate which is valid in the IHD limit. 

However, since the escape rate is simply given by the ratio of the total current through the saddle point to the 
number of particles in the metastable state, Langer's result for the escape rate can in fact be achieved by only 
computing the energy-Hessian eigenvalues near the saddle points and metastable states, from which one then infers 
the partition function Z s of the system restricted to the region around the saddle point where the energy- Hessian 
negative eigenvalue is (formally) 13] taken with absolute value, and the partition function Z m of the region around 
the metastable state. When computing these partition functions, one has to identify and take care of each Goldstone 
mode, that is a massless mode or zero-energy fluctuation associated with a continuous unbroken global symmetry. 
Finally, one computes the unique 01 negative eigenvalue k of the SSFPE corresponding to the unstable mode at 
the saddle point as the negative eigenvalue of the dynamic matrix M mn = —\ n (DMD T ) mn , where the A„'s are the 
eigenvalues of the Hessian at the saddle point and D is the transformation matrix from rji to ip n . 



4 



Consequently, Langer's final expression for the escape rate is rewritten in the following somewhat more practical 
form 

r = ^f, do) 

l-K Zj m 

where \k\ is the attempt frequency which contains the damping parameter a. 



III. RELAXATION RATE OF THE TWO-SPIN SYSTEM 



Now, we give the different relaxation rates for j > j c ,j < j c , and j ~ j c . 

Relaxation rate for j > j c : The quadratic expansion of the energy Q at the saddle point (0 reads, 

ffi 2) = 4 o) - |e + + f e, in) 

with zero eigenvalue for the £+ mode, that is the Goldstone mode associated with the rotation around the easy axis. 
ei°^ is given in Eq. The negative eigenvalue of the SSFPE, corresponding to the unstable mode, is k — —aj c . 

The partition function at the saddle point Z s = J dQidfl2S" l3Hs , where dQi = sin9id9idtpi, i = 1,2 is calculated by 
changing to the variables £±, C±, setting sin6*i ~ (see Eq. (0) in the integration measure, and finally computing 
the Gaussian integrals. Hence, 

^=2 7 r(^) 3/2 _=l=e-^, (12) 

\ I 3 ) jV 1 -Jc/j 

where we have formally replaced the negative eigenvalue (— j c ) by its absolute value. The partition function Z m at 
the metastable state © is computed by expanding the energy up to 2 n< ^ order, leading to 

Using Eq. 1)11) [) and inserting the symmetry factor of 2, the relaxation rate for j > j c reads 

T j>jc = 2aJf(l h) l + { y h),j x e-^ <0 \ AeW = (1 - hf . (14) 

For j ^ oo 1)14(1 tends to the Neel-Brown result, 

T j>jc -> 2«yj(l - h*){l - h) e-^\ (15) 

for the relaxation rate of one rigid pair of spins with a barrier height twice that of one spin. Note, however, that 
this convergence is very slow so that Tj > j c remains above the Neel-Brown result and only merges with the latter for 
j > 10. 

Relaxation rate for j < j c : Here we are faced with a semi-analytical case since the attempt frequencies \k £ \ (for 
e = ±) cannot be obtained in a closed form and are thus computed numerically. On the other hand, following the 
same procedure as for j > j c we obtain the relaxation rate for e = ±, 



3 ^R £ + R e _+jQ 6 {R+ + R e - 



where Ae^ is the barrier height given by the energy measured with respect to @ or Q for e = +, — respectively, 
and (seeEq. © et seq. for notation) P e = + b s h - j/2) 2 - (a e + h) 2 , Q e = b e h~j/2 + P £ ,R £ ± = -1 + X|(2X| + 
h), N + = (1 — h)(j + 1 — h), N~ = j c — j. The limit of the relaxation rate (|16|) when j — ► is just the Neel-Brown 
result for one spin. Indeed, the product of the last two factors in the prefactor tend to (1 — h)/«/2, the attempt 
frequencies tend to aj c = a(l — h 2 ), and the energy barriers Aef^ — > (1 — h) 2 /2. 
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FIG. 2: Left: Energy barrier as a function of j. Right: Relaxation rate (in logarithmic scale) of the two-spin system in the 
case of intermediate-to-high damping. The straight horizontal line is the Neel-Brown result with the barrier height taken twice 
that of a single spin. 

Note that in the present regime of j < j c the large value of anisotropy has not changed the temperature dependence 
of the individual relaxation rates, i.e., T^ <Jc , with e = ±. This is due to the fact that 1/yT appears in the prefactor 
each time there is a continuously degenerate class of saddle points 0, which is indeed the case for j > j c and j < j c 
with e = + and e = — . However, anisotropy do affect the temperature dependence of the relaxation rate of the 
two-spin system, since for j < j c there are two saddle points bringing each a factor 1/VT, see the first line in Eq. JSJ. 

The case of j ~ j c : When j approaches j c either from above or from below, more Hessian eigenvalues (in addition 
to A^+ ) vanish, rendering the saddle point rather flat and thus leading to a divergent relaxation rate. Indeed, for j ~ j c 
the relaxation rate (|14|l diverges, which clearly shows that Langer's approach which uses a quadratic approximation 
for the energy at the saddle point, e.g., Eq. Qllfl. fails in this case. The remedy is to push the energy expansion to 
the 6 th — order in the variable £_ (since A^_ = j — j c — ► as j — > j c ), i.e., 

5e s = e s -e^c^e + - 4 £ + ^£, (17) 

where c = (1 - j - 7/i 2 /4)/3 < 0, d = (j - 1 + 31/i 2 /16)/30 > 0. Then, the contribution ^2^/(3 of the mode 
£_ to the relaxation rate l|14|l must be replaced byf^° d£ > -e~^ Se , upon which the divergence of Tj>j c is cut off (see 
Fig. [21 (right)). Similarly, when j approaches j c from below, for e — — the eigenvalue A^- vanishes at j — j c , upon 
which the saddle point © (e = — ) merges with the state @ and thereby the partition functions Z s and Z m tend to 
infinity. However, as Z m oc l/(j c — j) increases much faster than Z s , and \k~\ — > 0, ^J < j c tends to zero as j — > j c . 
On the other hand, for e = +, both A^- and A^- vanish leading to a divergent relaxation rate, since now Z s diverges 
but Z m in Eq. I|13ll remains finite for the metastable state is well defined. Indeed, as j — > j Cl the relaxation rate 
Fj < j c goes over to the result in Eq. I|14l) upon making the change j <-> j c , and taking account of the symmetry factor. 

In this case, the divergence at the point j = j c cannot be cut off by expanding the energy beyond the 2 n( ^ order and 
^j<ic 1S s i m ply cu f onC & t t ne point where it joins Tj > j c taking account of Eq. I|17|l. 

IV. DISCUSSION AND CONCLUSION 

In Fig. |3 (right) we plot (In of) the relaxation rate of the two-spin system as defined in Eq. JSJ, under the condition 
of IHD, in which Langer's approach is valid, that the reduced barrier height /3Ae^ ^ 3> 1 and afiAe^ > 1 0- We 
also plot separately both relaxation rates for j < j c . We see that the relaxation rate of the two-spin system contains 
two unconnected branches corresponding to the two regimes, j < j c and j > j c , the bridging of which would require 
a more sophisticated approach. Fig. [21 (right) also shows that as j increases, but j <C j c , the relaxation rate Ejxj 
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decreases because the switching of the first spin is hindered by the (ferromagnetic) exchange coupling. While Tj<j is 
an increasing function of j with a faster rate, since now the exchange coupling works in favor of the switching of the 
second spin. This is also illustrated by the evolution of the energy barrier height in Fig. (left). As j approaches j c 
from below, the relaxation rate Tf<j c tends to Tj>j c because the respective saddle points merge at j — j c . Whereas 
^J<j § oes t° zero since the corresponding saddle point merges with the antiferromagnetic state that is no longer 
accessible to the system. For j > j c , as j increases the minimum ©, the metastable state © and the saddle point 
(JSJ) merge together, which means that the system is found in an "energy groove" along the direction 9\ = O2 because 
the eigenvalue corresponding to the mode 9\ — (9 2 becomes very large, and thus the escape rate decreases and 
eventually reaches the Neel-Brown value at large j. 

The present study has helped us understand the effect of exchange coupling on the relaxation rate of the two-spin 
system, and will be very useful for the generalization to multiple-spin small particles at least for small deviations from 
collinearity, where it has been shown |2( that the surface contribution to the macroscopic energy has a simple cubic 
anisotropy. However, this generalization can only be performed using numerical techniques. This is now attempted 
by the help of the ridge method ,101 f° r probing the energyscape and locating the saddle points, and by the (Onsager- 
Machlup) path integrals f° r determining the most probable paths connecting a metastable state to a more stable 
state. 
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